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LONG-TERM  GOALS 

1 .  We  seek  to  investigate  the  range  of  near-surface  mean  wind  directions  with  respect  to  the 
orientation  of  common  marine  atmospheric  boundary  layer  (MABL)  quasi-two  dimensional 
phenomena  seen  in  synthetic  aperture  radar  (SAR)  imagery. 

2.  We  seek  to  investigate  and  correct  for  the  impact  of  position  errors  in  numerical  model  synoptic 
scale  cyclones,  anticyclones,  fronts,  and  col  regions  on  the  wind  direction-dependent  retrieval  of  wind 
speed  from  SAR. 

OBJECTIVES 

1 .  We  will  investigate  the  range  of  near-surface  mean  wind  directions  with  respect  to  the  orientation 
of  common  MABL  quasi-two  dimensional  phenomena  seen  in  SAR  imagery  of  the  sea  surface.  We 
will  provide  empirically-derived  guidelines  on  how  to  distinguish  the  SAR-signature  of  one  feature 
from  another  (see  Figure  1  for  examples  of  quasi-two  dimensional  phenomena  seen  in  SAR  imagery  of 
the  sea  surface).  The  quasi-two  dimensional  MABL  phenomena  that  we  will  study  are  elongated 
cellular  convection,  buoyancy-driven  /  shear-organized  roll  vortices,  inflection  point-induced  roll 
vortices,  shear-driven  gravity  waves,  topographically-driven  gravity  waves,  and  dendritic  convection. 

2.  We  will  provide  the  impact  of  position  errors  in  numerical  model  synoptic  scale  cyclones, 
anticyclones,  fronts,  and  col  regions  on  the  wind  direction-dependent  retrieval  of  wind  speed  from 
SAR  via  CMOD-4.  See  Figure  2  for  an  example  of  such  an  error.  We  will  use  the  results  of  this  error 
analysis  to  develop  automated  and  semi-automated  detection  and  morphing  correction  algorithms  for 
displaced  synoptic  scale  features.  In  addition,  a  manual  editor  algorithm  will  be  developed  to  correct 
for  the  case  of  important  synoptic  features  missing  or  erroneously  included  in  the  numerical  model 
wind  direction  field.  These  algorithms  will  improve  the  accuracy  of  SAR  wind  speed  retrieval  in 
conjunction  with  numerical  model  wind  direction  fields. 

APPROACH 

1 .  The  basis  of  the  feature-identification  research  will  be  SAR  imagery  provided  to  us  by  Chris 
Wackerman  at  General  Dynamics  Advance  Information  Systems  (GDAIS).  GDAIS  has  a  collection  of 
SAR  imagery  from  which  linear  features  have  been  used  to  identify  corresponding  wind  directions. 
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Figure  1.  Radarsat-1  SAR  image  (left)  depicting  the  signature  of  roll  vortices.  The  300  m  pixel 
image  is  approximately  270  km  by  270  km.  The  image  was  acquired  at  C-band,  horizontal 
polarization,  off  the  north  east  coast  of  the  United  States  at  2242  UTC  on  6  March  1997.  The  top  of 
the  image  is  directed  towards  348° T.  (Provided  courtesy  of  JHUAPL,  ©  CSA)  ERS-1  SAR  image 
(right)  depicting  the  signature  atmospheric  gravity  waves  associated  highly  ageostrophic  flow  near  a 
front.  The  180  m  pixel  image  is  approximately  90  km  by  90  km.  The  image  was  acquired  at  C-band, 
vertical  polarization,  over  the  Caspian  Sea  at  0723  UTC  on  12  May  1996.  The  top  of  the  image  is 
directed  towards  012°  T.  (Provided  courtesy  of  Werner  Alpers  and  ESA,  ©ESA) 


Figure  2.  SAR  wind  speed  image  for  10  February  2002  generated  by  the  Johns  Hopkins  University 
Applied  Physics  Laboratory.  The  image  on  was  generated  using  NOGAPS  wind  directions.  The 
image  contain  an  example  of  an  erroneous  hour-glass  signature  (at  approximately  55°  N,  -145°) 
resulting  from  the  mis-positioning  of  the  cyclone  center  by  NOGAPS.  Arrows  denote  the 
corresponding  model  wind  speed  and  direction. 

We  plan  to  use  their  archive  to  assess  both  the  type  of  phenomena  depicted  in  a  given  SAR  scene  and 
the  relationship  between  its  orientation  and  the  actual  near-surface  mean  wind  direction.  Given  that 
hundreds  of  SAR  images  of  both  the  Northwest  Atlantic  Ocean  and  the  Northeast  Pacific  Ocean  are 
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available  to  us,  our  study  will  provide  a  robust  assessment  of  both  feature  orientation  with  respect  to 
near-surface  mean  wind  direction  and  of  how  to  discern  one  feature  from  another.  Using  this 
knowledge,  we  will  endeavor  to  develop  guidelines  to  distinguish  between  quasi-two  dimensional 
phenomena  based  on  their  structure  in  the  SAR  backscatter  field. 

2.  We  envision  three  methods  for  correcting  gross  errors  in  numerical  model  wind  direction  field  to 
obtain  a  revised  wind  direction  field  suitable  for  use  in  SAR-based  wind  speed  retrieval.  The  choice  of 
which  method  to  use  depends  primarily  upon  the  information  available  in  the  numerical  model  wind 
direction  field.  If  the  numerical  model  wind  direction  field  contains  all  of  the  synoptic  features 
observed  in  the  SAR-derived  wind  speed  field,  but  those  features  are  displaced  from  their  SAR 
positions,  then  an  automated  or  manually-guided  morphing  will  suffice  to  correct  the  position  errors. 

In  contrast,  if  important  synoptic  features  are  either  missing  or  erroneously  included  in  the  numerical 
model  wind  direction  field,  then  manual  editing  is  required. 

WORK  COMPLETED 

1 .  As  discussed  in  our  FY04  Progress  Report,  we  have  documented  the  SAR-signatures  of  synoptic- 
scale  cold  fronts,  warm  fronts,  occluded  fronts,  and  secluded  fronts,  including  their  mesoscale  and 
microscale  sub-structures.  In  doing  so,  we  ascertained  that  the  wind  direction  information  revealed  by 
the  SAR-signatures  of  mesoscale  and  microscale  eddies  is  of  great  importance  in  locating  and  typing 
synoptic  scale  fronts.  Thus,  there  is  a  direct  two-way  connection  between  this  front  research  and  both 
of  our  Objectives.  For  more  details  on  our  findings,  please  refer  to  Young  et  al.  (2005). 

This  front  research  has  resulted  in  a  follow-on  study  focused  on  distinguishing  between  oceanic  and 
atmospheric  fronts  in  SAR  imagery.  The  follow-on  study  is  being  conducted  by  Werner  Alpers, 

Jochen  Horstmann,  Andrei  Ivanov,  Leonid  Mitnik,  and  Todd  Sikora.  A  corresponding  journal  arcitle 
is  currently  be  prepared. 

In  conjunction  with  our  first  Objective,  we  have  discovered  a  new  form  of  quasi-two  dimensional 
boundary  layer  circulations  which  we  call  dendritic  convection.  These  branching  roll  circulations  form 
in  the  presence  of  cold  advection  over  tropical  oceans  when  the  wind  speed  is  less  than  that  required 
for  buoyantly-driven  roll  vortices.  The  orientation  of  the  resulting  convective  lines  provides 
information  on  both  the  surface  wind  direction  and  the  boundary  layer  shear  vector.  We  are  currently 
working  to  quantify  the  relationship  between  branching  angle  and  wind  speed  so  as  to  complete  the 
journal  article  that  is  currently  being  prepared. 

2.  We  have  completed  a  climatology  of  the  error  signatures  seen  in  SAR  wind  speed  imagery  as  a 
function  of  errors  numerical  model  fronts,  cyclones,  anticyclones,  and  col  regions.  Matlab  code  was 
developed  to  recognize  these  signatures  and  to  execute  manual  and  automated  single-control-point 
morphing  of  cyclone,  anticyclone,  and  col  region  positions  and  manual  multi-control-point  morphing 
of  frontal  positions.  In  addition,  Matlab  code  was  developed  to  enable  manual  editing  of  the  input 
wind  direction  field  for  the  case  of  synoptic  features  missing  from  the  numerical  model  wind  direction 
field.  A  corresponding  journal  article  is  currently  being  prepared. 

This  research  has  resulted  in  a  follow-on  study  focused  on  the  generation  of  a  SAR-derived  wind  speed 
risk  assessment  product.  Our  product  allows  one  to  assess  the  lower  and  upper  bounds  on  SAR-based 
wind  speed,  and  thus  the  wind-related  risk  to  maritime  operations.  The  risk  product  consists  of  a  pair 
of  wind  speed  images  bounding  the  wind  speed  range  consistent  with  the  observed  SAR  data.  We 
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envision  the  risk  product  to  be  most  useful  in  coastal  regions  where  the  large  concentration  of  maritime 
operations  requires  accurate,  high-resolution  wind  speed  data  and  when  uncertainty  in  the  a  priori 
knowledge  of  the  wind  direction  precludes  the  generation  of  accurate  SAR-based  wind  speed  fields. 

For  more  details  on  this  product,  please  refer  to  Sikora  et  al.  (in  press). 

RESULTS 

1.  We  have  documented  the  SAR-signatures  of  synoptic-scale  cold  fronts,  wann  fronts,  occluded 
fronts,  and  secluded  fronts,  including  their  wind  direction-related  mesoscale  and  microscale  sub¬ 
structures.  An  example  of  our  frontal  identification  findings  can  be  seen  in  Figure  3.  This  research  has 
been  formalized  in  Young  et  al.  (2005)  and  addresses  both  of  our  Objectives.  A  related  paper  on 
distinguishing  between  oceanic  and  atmospheric  fronts  is  being  prepared  by  Todd  Sikora  several  co¬ 
authors.  However,  much  work  remains  to  be  completed  in  connection  to  our  first  Objective.  The 
primary  obstacle  we  have  faced  is  securing  a  large  quantity  of  high-resolution  SAR  data  from  which 
the  signatures  of  linear  MABL  phenomena  can  be  discerned.  We  hope  GDAIS  SAR  collection  will  fill 
our  need. 

We  have  developed  a  space/time  climatology  of  dendritic  convection  and  a  conceptual  model 
explaining  its  structure  and  orientation  in  tenns  on  a  non-linear  interaction  between  buoyantly  driven 
quasi-two  dimensional  rolls.  This  conceptual  model  allows  us  to  determine  the  surface  wind  speed  and 
boundary  layer  shear  from  the  orientation  of  the  quasi-two  dimensional  structures.  Thus,  this  research 
relates  closely  to  our  first  Objective.  We  plan  to  complete  this  work  by  statistically  relating  the  surface 
wind  speed  to  branching  angle.  Our  research  is  being  formalized  by  a  journal  article  currently  in 
preparation. 

2.  We  have  developed  manual  and  automated  morphing  Matlab  algorithms  as  well  as  a  manual  wind 
direction  Matlab  editor  algorithm  in  fulfillment  of  our  second  objective.  We  have  used  our  Matlab 
products  to  demonstrate  that  substantial  improvements  in  the  SAR-derived  wind  speed  field  can  be 
obtained  by  correctly  editing  or  morphing  the  numerical  model  wind  direction  fields  before  application 
of  CMOD-4.  Our  research  is  being  formalized  in  a  journal  article  currently  in  preparation.  A  related 
follow-on  study  focusing  on  the  generation  of  a  SAR-derived  wind  speed  risk  assessment  product  has 
been  completed  and  is  fonnalized  in  Sikora  et  al.  (in  press). 

IMPACT/APPLICATIONS 

The  research  described  herein  fulfills  ONR  objectives  by  working  towards  the  integration  of  standard 
meteorological  theory  and  synthetic  aperture  radar  data  with  the  goal  of  providing  high-resolution 
remotely  sensed  estimates  of  near-surface  wind  field  in  in  situ  data-sparse  regions  over  the  ocean, 
including  the  coastal  zone.  For  example,  the  frontal  identification  and  risk  assessment  work  we’ve 
completed  thus  far  will  help  train  forecasters  to  use  SAR  to  its  full  potential  in  marine  synoptic  scale 
frontal  analysis.  The  morphing  software  we’ve  developed  will  allow  users  to  obtain  accurate  SAR- 
derived  wind  speeds  even  when  numerical  model  analyses  used  for  the  wind  direction  field  have 
position  errors  for  key  synoptic  features. 
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RELATED  PROJECTS 


George  Young  has  an  NSF-funded  SAR  project  that  supports  the  development  of  analysis  tools 
essential  to  our  ONR-funded  research.  The  Matlab  CMOD-4  code  was  developed  jointly  with  NSF- 
co-PI  Nathaniel  Winstead.  An  interactive  SAR  wind  field  analysis  program  with  a  graphical  user 
interface  was  developed  by  George  Young  as  part  of  this  project  and  will  serve  as  the  host  application 
for  the  automated  and  interactive  wind-direction  field  morphing  routines  developed  under  ONR 
funding. 
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Figure  3.  Radarsat-1  SAR  image  depicting  the  the  intersection  of  a  cold  front,  warm  front,  and 
occluded  front.  The  600  m  pixel  image  is  approximately  800  km  by  450  km.  The  image  was 
acquired  at  C-band,  horizontal  polarization,  over  the  Bering  Sea  at  0500  UTC  on  31  October  1999. 
The  top  of  the  SAR  image  is  directed  towards  348°  T.  (Provided  courtesy  of  JHUAPL,  ©  CSA) 


Todd  Sikora  is  collaborating  with  Werner  Alpers,  Jochen  Horstmann,  Andrei  Ivanov,  and  Leonid 
Mitnik,  on  a  study  focused  on  distinguishing  between  oceanic  and  atmospheric  fronts  in  SAR  imagery. 
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